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INTRODUCTION 


Results of previous NAS/l support showed that multi spectral scanning tech- 
nology was capable of measuring the concentration of chlorophyll In benthic 
microalgae (mainly diatoms) on an estuarine mudflat (c.f. Final Report of NASA. 
Grant No. NSG-1334, 1978). Modest support was requested In FY 79 to conclude 
work Initiated in NSG-1334. Specifically, we needed additional time to 1. 
analyze further the large volume of data that we had accumuJated on magnetic 
tape; and 2. to communicate the results of our research to the scientific 
community. 

DATA ANALYSES 

The data we collected In July 77 was Initially very puzzling, as the 
correlation coefficients between the ground truth and the scanner derived 
data were disturbingly low (Table 1). We determined the cause of the poor 
correlations through an in-depth analysis of the data. 

An intertidal estuarine mudflat Is not uniformly flat. Rather, such 
muflats are alternately covered and uncovered with flowing tidal water twice 
each day in southeastern estuaries. The movement of the ebbing tide tends 
to form small, anastlmoslng rivulets that get larger as they approach the center 
of the mudflat. These rivulets lie between slightly elevated ridges. Col- 
lectively, the ridges and rivulets form a topographically complex surface In 
contrast to tlie rclntivclv uniformly flat surface of the remainder of the 
onudflat. It was obvious from the graphically represented data that some of 
the data was clustered and other was scattered (Fig. 1). Thus, we separated 
. the data derived from the rivulet areas from that from the flat areas and re- 
calculated the correlation coefficients. This time the ground truth data from 


Table 1. Correlation coefficients of all ground 
truth chlorophyll measurements and coincident 
oignals levels of the multispectra 1 scanner in 
July 1977. 

SPECTRAL 

CHANNEL 

CORRELATION 

COEFFICIENT 

BLUE 

-0.20 

GREEN 

-0.125 

RED 

-0.155 

IR 1 

0. 13 

IR 2 

0.33 

IR 3 

0.36 


Che flat areas correlated well with chat of the blue channel of the scanner, 
while the data from Che rivulet areas was only slightly better correlated 
(Table 2). Graphic representations of the data make the differences stand 
.out clearer (Figs. 1, 2). We conclude the obvious necessity to partition 
future data into subsets prior to computing the concentration of chlorophyll 
over the whole mud flat. 

DISSEMIN/iTION OF RESULTS 

I attended and presented the results and conclusions of this research 
at Che Fifth International Symposium of Living and Fossil Diatoms in Antwerp, 
Belgium (Appendix 1). The paper was well received and many fruitful discus- 
sions with colleagues followed. The paper was accepted for publication by 
the editor of Che symposium's proceedings, but I declined the invitiation. 

I judged that Che paper would be read by a larger and wider audience and be 
published sooner, if sent to a remote sensing journal. Following the symposium 
I visited the laboratories of Drs. Wim Admiraal and F. Colijn of the University 
of GrBningen, Netherlands, and with Dr. Klaus Wegmann of the University of 
TUbingen, West Germany. All are engaged in remote sensing research on benthic 
diatoms. These were rare and valuable opportunities to discuss in detail the 
import of our research and the future of multlspecCral scanning technology in 
estuarine and coastal ecological research. 

The initial results of our research Iiave been submitted to the journal 
Remote Sensing of the Environment and are presently undergoing peer review • 
(Appendix 2). 

CONCLUSIONS 

We have determined that there is a stastically significant relationship 
between the ground truth measurement of chlorophyll of benthic microalgae and 
•radiance levels in the blue spectral channel' of a tower-mounted multispectra 1 


Table 2. Correlation coefficients of subsets of all .ground 
truth chlorophyll measurements and coincident signal levels 
of the multispectra 1 scanner in July 1977. 


SPECTRAL 

CORRELATION COEFFICIENTS 

CHANNEL 

RIVULET AREA 

FLAT AREA 

BLUE 

- 0.33 

-0.82 

GREEN 

- 0.45 

-0.37 

RED 

-0.46 

-0.24 

IR 1 

-0.02 

0.25 

IR2 

-0.18 

0.51 

IR3 

0.18 

0.36 


GROUND TRUTH 



Figure 1. Scattergram plot of all ground truth chlorophyll' 
measurements and coincident signal levels of the multispectra 1 
scanner in July 1977. 
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Figure 2. Replot of the scattergram seen In Fig. 1 but Indi- 
cating the significant difference between the subsets. 



scanner, provided that the data is first partitioned Into similar subsets. We 
are encouraged by these results, have reported on them at scientific meetings 
and have submitted them for publication In the open scientific literature. We 
anticipate that continued financial support by N^SA will allow us to experi- 
ment with aircraft remote sensing of the benthic microalga 1 connmnlty, an 
ecologically significant coniminlty In estuarine ecosystems. 


APPENDIX 1. Abstract of a paper read at the Fifth International Symposium on 
Living and Fossil Diatoms, Antwerp, Belgium, September, 1978 
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Abstract: 

A remote sensing instrument, similar in design to the multispectra 1 
scanner on the NASA VIKING lander on Mars, was mounted on a 50 foot tower 
overlooking North Inlet Estuary, South Carolina. The instrument was pro- 
grammed to take multispectra 1 imagery data along a 90° horizontal frame, 
while measuring relative radiance in six spectral bands ranging from 400- 
1100 nm and had a ground resolution of 2-5 cm. Imagery measurements were 
encoded in digital form on magnetic tape and were stored, decoded, and ma- 
nipulated by computer. Correlation coefficients from data on scans of ti- 
dally exposed mud were highest in the blue and were negative, possibly in- 
dicating the absorhtion of sunlight by the chlorophyll containing benthic 
microflora (mainly diatoms). Concurrent, quantitative "ground truth" mea- 
• surements of extracted chlorophyll a from cores were made to calibrate the 
digital data as recorded by the scanner. The data from the two widely dif- 
' ferent techniques had correlation coefficients between 0.81 and 0.91 on 
three separate sampling periods. Seasonal patterns of chlorophyll concen- 
tration on the nmdflat followed a predictable pattern, with a winter low 
and a spring high. The scanner has provided encouraging results and pro- 
,mises to be a useful tool in sampling the biomass of the benthic micro- 
flora over large intertidal estuarine areas. 


APPENDIX 2. Copy of tnanuscrlpt sent to the journal. Remote Sensing 
of the Environment. ' 
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ABSTRACT 


A remote sensing instrument, vas mounted on a 50 foo> tower overlooking 
North Inlet Estuary, South Carolina. The instrument was programmed to take 
multispectral imagery data along a 90° horizontal frame in six spectral bands 
ranging from UOO-1050 nm and had a ground resolxition of about 3 cm. Imagery 
measurements were encoded in digital form on magnetic tape and were stored, 
decoded, and manipulated by computer. Correlation coefficients were 
calculated on imagery data and chlorophyll & concentrations derived from 
grovind truth data. The most significant correlation occurred in the blue 
spectral band with numerical values ranging from -0.8l to -0.88 for three 
separate sampling periods. Mean values of chlorophyll a for a larger section 
of mudflat were estimated using regression equations. The scanner has 
provided encoiiraging results and promises to be a usefvil tool in sampling 
the biomass of intertidal benthic microalgae. 


INTRODUCTION 


Benthic microalgae are significant components of estuairine food webs 
(Pomeroy, 1959 i Marshall, Oviatt and Skauen, 1972; Ranwell, 1972), and 
probably play an important role in stabilizing estuarine sediments (e.g. , 
Holland, Zingmark and Dean, 197^)- Yet most of the literature on primary 
producers in estuaries has concerned the vascvilar plants and the phytoplanktonic 
commiinities. Consequently, ovir present understanding of benthic microalgal 
biomaas and productivity is based on a relatively small number of data points. 

The biomass of the benthic microflora is usually estimated by measvuring 
chlorophyll a in discrete cores of svirface sediments; (Grjintved, I960, I 962 , 
1965 ; Gargas, 1970; Cadee and Hegeman 197*0- The number of such samples 
that can be taken in any area is limited by time and financial constraints. 

Yet, because the benthic microflora does not form a homogeneous layer, a 
large nmber of replicate samples must be measvured to provide a ’’representative 
sainple" of an area. This is especially true when one considers that the 
composition and size of the sediments influences to a large degree the 
species composition and biomass of the benthic microflora (Sullivan, 1975; 
Ainspoker and MacIntyre, 1977; DeFelice and Lynts, 1978) and that the 
composition of estuarine sediments is variable depending on bottom topography 
and the velocities, current flow and wave patterns of the overlying water 
which influence the deposition, mixing, and erosion of sediments. 

More extensive sampling of benthic microalgal communities is necessaiy 
to understand the differences in biomass on the various sediment types than 


is financially feasible, especially in large estuarine areas. Frequent 
personal observations in North Inlet Estuary, South Carolina, at low tide 
have demonstrated that exposed benthic microalgae are often visible as 
gold to brownish patches or films to the unaided eye. These observations 
illustrate the potential of tising a remote sensing technique for measuring 
the chlorophyll biomass of this community. Such methods have the potential 
of sxirveying large areas in a short time at a reasonable cost. . 

Remote sensing of the primary converters of solar energy in estuarine 
and marine food webs have heretofore focused on measuring dissolved and 
particulate matter in the water, and mapping coastal vegetation and other 
terrestrial features (Clarke, Ewing and Lorenzen, 1970; Stoertz, Hemphill 
and Markle, 1970; Carter and Schubert, 197^)* These studies demonstrated 
the feasibility and usefulness of obtaining reliable, duplicable measurements 
that were largely independent of the sample size, since such methods scanned 
a wide surface area. 

We report here an investigation designed to determine the feasibility 
of using existing multispectral technology for detecting and recording high 
resolution, quantitative information on the biomass of benthic microalgae 
on an intertidal, estuarine raudflat in South Carolina. The primary purpose 
of this study was to determine whether a statistically significant relation- 
ship could be demonstrated between relative spectral radiance measured by a 
tower mounted, multispectral scanning instrument and the concentration of 
chlorophyll a measured in discrete ground truth samples . This approach was 
chosen primarily for reasons of high spatial resolution and low cost. The 
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grovind resolutions available from existing aircraft scanners ranges from 
about T to TO meters while the minimum ground resolution of the ctirrent 
LAIIDSAT multispectral scanner is about 80 meters. Current techniques of 
in-situ sampling use core samples a few centimeters in diameter. In view 
of the complete absence of any background research in remote sensing of 
benthic microalgae we considered an exact match in scale between image 
resolution and ground truth samples was essential. The tower-mounted 
scanner provided this match with a spatial resolution of approximately 3 
centimeters at the sample area. Likewise, economics was a prime factor in 
choosing a tower platform since aircraft, depending on size cost typically 
$150 to $1,000 per fli^t hoxir. The cost of operating the tower mounted 
scanner is negligible. While a tower mounted scanner sufficed for this 
feasibility study, an aircraft scanner would be required for research 
requiring large area coverage. 

The study site was a 10 x 200 m intertidal mudflat emptying into 

Clambank Creek, North Inlet Estuary, South Carolina (Figure l) . North 

2 

Inlet Estuary is a 30 Km , hi^ salinity, tidal, salt-marsh estuary 
dominated by the marine halophyte Spartina altemiflora . 

SCANNER SYSTEti DESCRIPTION 

Present day multispectral scanners unlike film cameras take image data 
as a sequence of image elements which form a line. Then a sequence of lines 
of image are taken to form a frame in each spectral band of interest. Tlie 
basic physical quantity measured by any multispectral scanner is spectral 
radiance (power per unit area - solid angle - wavelength interval) either 
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reflected or emitted by a scene. A sennner responds with an electrical sicnol 
proportional to the spectral radiance of an image clement defined by the 
instmitaneous field of view (IFOV) (See Figure 2). This ITOV is set by the 
size of a photodetector and its distmice from the imaging objective lens. 

For far field imaging this distimce is appi'oximately the focal length of 
the lens. 

The signal current due to the spectral radiance of an image element is 

1 = kA!! / nCA)t(X]R(X)c 1X (l) 


for a glvt'n spectral band AX = - X^ where k is a proportionality 

constant representing primarily elccti'onic gains and/or calibration factors, 
A is the aperture area of the imaging objectives, fi, the solid angle of 
the IFOV, N(X), thd spectral radiance of the scene element. t(X), the 
spectra! transmissivity of optics R(X), the spectral responsivity of the 
photodetector, and dX the diffcrontitil wavelength. 

For far field imaging witii a circular photodetector 




Cf.l.)‘ 


( 2 ) 


wliere d is the photodetector diiunoter and F.I,. , the lens focal length. 
'Ihe current signals are then usually converted to voltages and digitized 
for telemetry or direct recording on mag.netLc tape. 

Ibe i>articular multi spectral scanner used in this investigation was 
designed ?md built at NASA Iiangley Research Center except for the array of 
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photodetectors. A detector array which was a backup flight unit for the 
Viking Lander Cameras was used because it was readily available and had 
known spectral responsivities which provide the desired measurement 
capability. A list of scanner characteristics is given in Table 1. One 
item not listed is cost. It is of interest that the cost of this instrument 
would be lower, if the Viking detector array is excluded. This array was only 
used as a matter of convenience and could have been replaced by a more 
conventional detector array with interference filters at a cost of less 
than about $1,000. The resulting overall cost would have been under 
$ 11 , 000 . 

Since the scanner was designed to be operated by unskilled technr.cians, 
most of the scanner operation was automated to such an extent that the 
operator could take frames of multispectral image data by manipulating a few 
switches. The automatic control functions are scanning of lines, switching 
of detector elements to provide line sequential color and multispectral 
near infrared images, and data recording. All azimuth positioning to 
start a frame is perfoi*med automatically. The operator is required to 
select the extent of the azimuth frame , to turn the scanner controls and 
data storage device on and off, and to insert and replace the magnetic 
tape cassettes. 

The cassettes are the datci storage medium for the digitized multi- 
spectral data. This mediiun was selected for compactness, light weight, 
ease of operation, and ease of shipment of the cassettes. Ihis later 
reason was important since image data analysis was performed in Virginia. 

The cost of the data storage unit was $2500. 
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EXPERIMENTAL PROCEDURE 


All experiments vere run on completely clear or uniformly overcast days 
to eliminate signal variations due to moAnng clouds. Dates of sampling were 
chosen when a minus low tiae was predicted between 10 a.m. and 2 p.m. to 
provide high sun angle conditions. 

The experimental procediare began by manning the scanner through its 

programmed cycle after the mudflat had been exposed to dry for at least 

one hour. Following this initial run, one or two linear series of at 

2 

least eighteen 5«^3 cm x 1 cm cores were taken for chlorophyll analysis. 

The depression left by the removed mud was covered by a 200 ml white 
styrofoam cup to mark the specific image elements where ground truth 
measurements were taken. This was followed by running the scanner through 
its cycle a second time to record the specific sites of ground truth samples 
on magnetic tape. Later it was possible to calibrate the scanner by 
comparing the radiance values as measured by the scanner at the ground 
truth sites with the laboratory values obtained for chlorophyll from 
these sites. Then it was possible to estimate the distribution of chlorophyll 
over the whole mudflab transect. 

Chlorophyll was extracted from well mixed, moist cylindrical core 
samples (2.5^ cm dia. x 1 cm deep) in 25 ml 90 percent Acetone in the dark 
at for four hours. Subsamples of 0.05 to 0.10 ml from each extract was 
diluted to 10 ml with 90 percent Acetone and the resulting mixture analyzed 
^ for chlorophyll a in a Turner JSTR Fluororaeter (Strickland and Parson, 19YT)* 

mu 
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DATA ANALYSIS 


A brief discussion of radiometry is given to provide an understanding 
of the variables which influsnce the multispectral imagery data of the 
mudflat. In addition, this discussion centers on the inability to 
determine benthic microalgae chlorophyll concentrations directD.y by 
radiometric analysis. The statistical approach which was used instead 
is detailed. 

Radiometric Considerations 

The spectral radiance sensed by a multispectral scanner for a mudflat 
is a combination of diffusely reflected direct sunlif^t and specularly 
reflected skylight. The reflectance properties of the mudflat are also 
complicated by the algae film-miud substrate conibination . The reflectance 
o'.’ •fne mud is a function of surface structure, topography, moisture content 
and soil composition. Even if all these variables could be taken into 
account the chief hindrance to radiometric analysis would still be the 
unknown relationship between the optical properties of the algap film 
and the chlorophyll concentration of the film. Since radiometric analysis 
poses seemingly intractable problems, a straightforward statistical approach 
is used. Two aspects of radiometry however are still pertinent to a 
statistical analysis of the multispectral data, since they deal with sources 
of extraneous variability. In this case the two factors are terrain 
features and view angle effects. The main terrain feature is a "rivulet 
pattern" confined mainly to the center portion of the mud flat. The view 
angle is a variation in signal with view angle which is due to the specular 



corroonent of mu'! flat reflectance being a fairly strong function of view 
angle for the vertical view angle range of the multispectral scanner. The 
two factors are eliminated as sources of extraneous variability in the 
following way: 

(1) samples occvurring in "rivulet patterns" are excluded from the 
analysis, and 

(2) all samples are taken in horizouteO. lines of nearly constant view 
angle. 

Correlation 

A correlation analysis weis performed on chlorophyll concentrations 
and the multispectral signals in the six visible and near infrared 
spectral bands for image elements where chlorophyll concentrations were 
measured. This analysis takes the usual form. 


m 
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correlation analysis was based upon the appearance of a roughly linear 
relationship in the data itself. This type of relationship cannot be 


expected in general. 

Since the statistical distribution of r is not normal, a test 
of significance of r is on Fisher's transformation (Freund, 19o2) 
of r 



which does have an approximately normal distribution. The Z statistic 
can then be used in a one tailed test to deteimine the miiMu.ura value of 
that allows the hypothesis |r| £ |r^| » to be rejected by 



^ - k 


N is the number 
This test can be 
confidence ?.evel 


of samples and k the number of estimated parameters, 
applied for any confidence level but for 95 percent 
z > 1.61|. 


Regression 


The goal of using remote sensing data to enlarge the number of samples 
of some variable compared to direct sampling can be performed with some 
confidence if a significant correlation is found between scene spectral 
radiance and the independent variable of interest. In this case an 
estimate of a regression equation is developed as 


c = ao + 


from the small sample by (Miller and Freund, 1977) where 
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Using equation (7) a population mean of chlorophyll concentration can 
be estimated from camera signals. For large populations, the dominant ^ 

\ 

error in absence of extraneous variabilities in the scene is in the 
estimates of a^ and a^ since these errors are systematic biases in 
the population mean estimate of chlorophyll concentration. The bounds 
on a^ and a^ estimates can be determined for a given confidence level 
(Miller and Freund, 1977) as 
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where t is the t statistic for 1 - a confidence level and 
a/2 
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The bounds which are calculated from (lO) and (ll) are then used to 
determine two lines which bound the regression line for a given confidence 
level. 


RESULTS AND DISCUSSION 


Multi spectral image data and ground truth samples were taken on 
three dates May 2, Jtily 26, and September 23, 1977* Data taking oppor- 
tunities were limited by the number of days when low tides occurred at 
appropriate times and by weather. 

Correlation Coefficient 

The correlation coefficients for the three dates are given in Table 2 
together with the number of samples. For the July 26 data some samples 
were taken in the "rivulet" pattern of the mud flat. These samples were 
eliminated from the correlation analysis. The September 23, 1977, data 
contained some areas with direc Fun specular reflections. These were 
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also eliminated from the correlation analysis hy a threshold criterion. The 
results show a consistent high negative correlation in the blue spectral 
band. Scattergrams are presented in Figure h for the blue spectral band 
for all three dates. The high negative correlations in the green and red 
spectral bands for the September data cannot be explained at this time. 
Possible explanations are that the dominant species of microalgae may have 
shifted to some new species with different spectral properties or possibly 
some strictly seasonal change has occ\irred without a change in species. 

Hypothesis Testing 

The significance test used here is a one tailed test to determine 
how much lower the population correlation coefficient could be than the 
sample correlation coefficient and still maintain a 95 percent confidence 
in rejecting the hypothesis that the population correlation coefficient 
was lower than that value. Values for Table III are given for which 
z = 1.64 indicating the magnitude of the population correlation 
coefficient was equal or greater than |r^| with a 95 percent confidence 
level. This result confirms that a significant negative correlation was 
fovind in all three sets of data. 

Linear Regression 

Since the correlation coefficient for the blue is consistently and 
significantly below zero, a linear regression relationship can be derived 
with some confidence. This regression relationship can be expected to 
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have limited application for the mud flat. It cannot he considered to 
describe view angle effects or rivulet pattern variations. As a result 
blue image elements signals can be used to esti’ e chlorophyll concen- 
tration of benthic microalgae only over a transect region of the mud 
flat for which samples were taken, and within this s\ibset of data, the 
rivulet pattern pixels were excluded. With these restrictions a 
regression relationship was used to estimate the mean chlorophyll concen- 
tration for two or three thousand image elements based on the 8 to 20 
core samples- This application of each regression equation yielded 
estimates 'of transect mean values of chlorophyll concentration which are 
given in Table U together with core samples means presented for comparison. 

Discussion 

<\ 

The concentration of chlorophyll pi as measured from discrete ground 
truth samples varied widely from sample to sample, indicating a 
patchiness of distribution of the benthic microflora. May samples 

A _p _ 

varied from U.09 - 7*79 mg Chid • m (x = 5*02), July samples from 1.8i* - 
(y _o _ 

6.32 mg Chid • m (x = U.07), September sr-mples from 0.86 - 7*50 
^ “2 ■" 

mgClild • m (x = 3.5^)- A steady decline of chlorophyll concentration 
was apparent throng the course of this study as seen from the ground 
truth data. 

Using the regression equations, we estimated the mean chlorophyll 

concentration over a large block of mudflat. Mean values thus calculated 

—2 ' —2 

were U .65 mg Chid • ra CMay), 5.65 mgChld • m (July) and 3-^0 mgChla • m 
(Sept.), using from 1900 — 2700 pixel elements. A peak in chlorophyll 


was seen in the Summer in the data recorded by the scanner. By placing 
the means on the regression curves we could estimate 95 percent confi-- 
dence intervals in the remote sensing data (Figure U, Table 2). 

Previous studies have indicated that the chlorophyll biomass of 
the benthic microflora fluctuates seasonally, with lowest values seen 
in the cold winter months and hi^er values seen during the warmer 
months (e.g. Marshall et al. , 1971)* However, sharp seasonal peaks, so 
well known for phytoplankton, are not as well defined (Cadee and Hegeman, 
197^)* Our data were taken at too infrequent intervals to show any distinc- 
tive seasonal pattern. However, there are differences between the 
pattern for the mean values for chlorophyll in the two methods used in 
this study. Highest values for chlorophyll were measured in the cojre 
samples in May, while the highest values sensed by the scanner were in 
July (Table 2). The variations in chlorophyll concentration between each 
of the discrete core samples was very great, and the discrepancy between 
the means for the two methods may be due to the relatively small sample 
set represented by the ground truth^e^. 

C ONCLUSIONS 

A statistically significant linear inverse relationship was found 
between chlorophyll concentration of benthic microalgae and radiance 
levels in a blue spectral band as measured by a tower mounted multispectral 
scanner. This result was found for three sets of data triken at different 
times during the growing season. One set of data exhibited high negative 
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correlation coefficients in the green and red spectral bands as veil. In all 
data, two major sources of extraneous variability had to be excluded 
from the correlation analysis. These were view angle effects and a 
’’rivulet" pattern in the center of the intertidal mud flat. These 
variables were eliminated by experiment planning and data selection. 

This statistical analysis of course cannot determine whether the 
relationship is a direct or an indirect one. The relationship, if direct, 
could be due to scattering and/or absorption by the algal film or, if 
indirect, possibly is a measure of moisture variations across the mud 
flat. 

We are encouraged by these results that aircraft remote sensing of 
benthic microalgae may be possible. Anticipated problems which were not 
encoxantered in the tower measurement are scale mismatch between ground 
truth and imagery picture elements, atmospheric path radiance effects, 
and additional extraneous svtrface variables such as sediment type varia- 
tions all of which may complicate an aircraft remote sensing measurement 
of benthic microalgae. 
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Comparison of ground truth neasurements of chlorophyll 
and estimates of chlorophyll as sensed by a multispec- 
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fiijiire Scanner schematic. 








Figure 3.- Intertidal mudllat showing sample area outlined 

in white. This image was reconstructed from tower 
scanner data. The two sets of non-inage lines occur 
during track switching by the cassette recorder. 







